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Abstract 

This investigation employs a square lattice and Glauber dynamics methodology 
to probe the effects of diluteness in the crystal-field interaction in a Blume-Capel 
Ising system under an oscillating magnetic-field. Fourteen different phase diagrams 
have been observed in temperature-magnetic-field space as the concentration of the 
crystal-field interactions is varied. Besides, a comparison is given with the results 
of the pure spin-1 Ising systems. 
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1 Introduction 



Ising spin-1 systems, with density as an added degree of freedom, have 
been utilizied to investigate a diverse range of systems: materials with mobile 
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defects, structural glasses [T] , the superfluid transition in He 3 — He 4 mixtures 
[2], frustrated Ising lattice gas systems [3], binary fluids, binary alloys and 
frustrated percolation |3]. Recently, Blume-Capel (BC) [5] and Blume- Emery- 
Griffiths (BEG) models have been used to probe the phenomenon of inverse 
melting, a phenomena observed in diverse class of systems such as colloids, 
polymers, miscelles, etc [5]. 

On the other hand, it is well known that the effects upon criticality and 
resulting phase diagrams, due to underlying competing interactions in various 
spin-1 Ising systems can be complicated. Since then, many previous investiga- 
tions have been focused on various novel types of competing interactions have 
been the focus of previous studies using the BEG model in conjunction with 
renormalization-group [7|S|9|TU] and/or mean-field methodologies [TTfl2|ll3] . 
Besides, effects of disorder on magnetic systems have been systematically stud- 
ied, not only for theoretical interests but also for the identifications with ex- 
perimental realizations [H|15II16] . It has been shown by renormalization group 
arguments that first-order transitions are replaced by continuous transition, 
consequently tricritical points and critical end points are depressed in temper- 
ature, and a finite amount of disorder will suppress them |17j . 

An special magnetic system with disorder is spin-1 Ising model with quenched 
diluted single ion anisotropy is used to model phase separations of superfluid- 
ity for helium mixtures in aerogel [THfTU] . Due to this fact, various researchers 
have been motivated to study the effect of the crystal field disorder on the 
multicritical phase diagram of BC model via effective field theory |20j and 
mean field approach [2T1T22] .cluster variation method |23j, as well as by intro- 
ducing an external random field [23]. Whereas, Branco et al. considered the 
effects of random crystal fields using real-space RG [8,9] and mean- field ap- 
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proximations [9~f25] for both BC and BEG model Hamiltonians, respectively. 
Recently, Snowman has employed a hierarchical lattice and renormalization- 
group methodology to probe the effects of competing crystal-field interactions 
in a BC model [26J . Finally, Salmon and Tapia have studied the multicriti- 
cal behavior of the BC model with infinite-range interactions by introducing 
quenched disorder in the crystal field A i; which is represented by a superpo- 
sition of two Gaussian distributions [27]. 

While the equilibrium properties of the BC model with random single ion 
anisotropy have been studied extensively, as far as we know, the kinetic aspects 
of the model have not been investigated via Glauber dynamics. Therefore, the 
purpose of the present paper is, to present a study of the kinetics of the spin-1 
BC model with a quenched two valued random crystal field in the presence of 
a time-dependent oscillating external magnetic field. We make use of Glauber- 
type stochastic dynamics to represent the time evolution of the system [28] . 
More precisely, we have obtained the dynamic phase transition (DPT) points 
and presented phase diagrams in constant crystal field and the reduced mag- 
netic field amplitude versus reduced temperature plane for various values of 
the crystal field concentration. This type of calculation for pure BC model, was 
first performed by Buendia and Machado [29J. They have presented only two 
phase diagrams in the temperature-magnetic field plane for the pure spin-1 
BC model. Later, Keskin et. al. have shown that one of the two phase dia- 
grams in Ref (29] was incomplete; i.e., they had missed a very important part 
of the phase diagram due to the reason that they did not make the calcula- 
tions for higher values of the amplitudes of the external oscillating magnetic 
field |30j. Keskin et. al. presented the phase diagrams in the reduced magnetic 
field amplitude (h) and reduced temperature (T) plane and calculated five dis- 
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tinct phase diagram topologies. Recently, Elyadari et.al. has investigated the 
kinetic Blume Capel Model with a random crystal field distributed according 
to the following law: P(Aj) = p5(A { - A(l + a)) + (1 - p)<J(Aj - A(l - a)) . 
This kind of random crystal field has been introduced to study the critical 
behavior of 3 He — 4 He mixtures in random media (aerogel) modeled by the 
spin-1 Blume-Capel model [31] . In their model, the negative crystal field value 
corresponds to the field at the pore-grain interface and the positive one is a 
bulk field that controls the concentration of 3 He atoms. In other words, in 
this kind of randomness the crystal field value is finite for each site about its 
amplitude takes one of the values A(l + a) or A(l — a) with equal and fixed 
probabilities (p = 1/2). While in our current investigation, we have focused 
on another kind of randomness in which one can see the effect of vacancy 
of the single-ion anisotropy on the dynamical phase diagrams of the kinetic 
Blume-Capel model. Previous equilibrium studies has revealed that this kind 
of randomness leads to a rich variety of phase diagrams with type being ac- 
cording to the concentration p of active local crystal fields. |20)21j . With this 
motivation we have performed numerical calculations for various values of the 
crystal field concentrations in order to observe the effect of the quenched va- 
cancy in the crystal field on the five different kinetic phase diagram topologies 
found by Keskin and co-workers [30] . 

Meanwhile, it is worthwhile to stress that the DPT was first found in the 
study of the kinetic Ising system in an oscillating field [32J, and it was fol- 
lowed by Monte Carlo simulation researches of kinetic Ising models [33~f34"] . 
Further, Tutu and Fujiwara [35] represented a systematic method for obtain- 
ing the phase diagrams in DPTs, and constructed a general theory of DPTs 
near the transition point based on a mean-field description, such as Landau's 
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general treatment of the equilibrium phase transitions. DPT may also have 
been observed experimentally in ultrathin Co films on Cu(OOl) [36] by means 
of the surface magneto-optic Kerr effect and in ferroic systems (ferromagnets, 
ferroelectrics and ferroelastics) with pinned domain walls [39] and ultra-thin 
[Co/Pt] 3 multilayer [37]. In addition, reviews of earlier research on the DPT 
and related phenomena are found in Ref [34J . 

The paper is organized as follows: In Sec. 2, we discuss the kinetic BC model 
with single ion isotropy briefly. Moreover, the derivation of the mean-field dy- 
namic equations of motion is given by using a master equation formalism in 
the presence of an oscillating external magnetic field is also given in Sec. 2. 
In Sec. 3, the DPT points are calculated and the phase diagrams presented. 
Finally, Sec. 4 represents the summary and conclusions. 



2 DYNAMICAL EQUATIONS FOR THE MEAN VALUES 

The generalization of the kinetic BC model for a quenched random crystal 
field is given by the Hamiltonian, 

H = -Jj2 S i S i ~ E A i s ? " , (1) 

(ij) {i} « 

where the spin variables Si = 0, ±1 on a square lattice. The first and the 
second sums are over nearest neighbor pairs. The exchange interaction with 
strength J > is responsible for the ferromagnetic ordering, while the random 
single ion anisotropy Aj is given by the following joint probability density: 

P{Ai) = p6(Ai - A) + (1 - p)S(A i ) . (2) 
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Finally, if is a time-dependent external oscillating magnetic field and given 

by, 



H{t) = H cos(ut) , 



(3) 



here Ho and u) = Inv denote the amplitude and the angular frequency of the 
oscillating field respectively. 

When we put this system in contact with a heat reservoir at temperature 
T, the spin variables Si can be considered as stochastic functions of time. The 
system evolves according to a Glauber-type stochastic process at a rate of - 
transitions per unit time. More precisely, we will follow the heat-bath prescrip- 
tion [38J: the new value of the spin variable at site i (S inew ) is determined by 
testing all its possible states in the heat-bath of its (fixed) neighbors (here 
four on a square lattice): 



where ft = and r defines a time scale (characteristic mean time interval for 
one spin flip), and 



give the changes in the energy of the system in the case of flipping of the i 
spin in the lattice. If we define P(Si, S 2 , S/v! t) as the probability that the 
system has the configuration {Si, S2, S^}, at time t. Making use of master 
equation formalism [28J, one can write the time derivative of P(Si, S2, Sjv; t) 




1 exp{-f3AE(Sj i d -» Sj new )} 
TT,exp{-f3AE(S io id -> S inew )} 



(4) 
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—P(Si, Sn; t) — —Y,i Y,Si old ^s inew Wi(Si id — > Si new) P (Si, Si old, ...Sn; t) 

+ J2i J^Si i d ^Si new w i(Si new Si ld)P(Sl, Si new, ...Sn', t) 

The detailed balance condition reads, 

(7) 



U!i(Si old Si new) _ P(Sl, S2, Si new , ■ Sn) 
Wi(Si 

new ^ Si old) P(Sl, S2, Si old, •••Sn) 



In addition, substituting the possible values of Si new and Si Q id , one ob- 
tains: 

Wi(l -> 0) = Wi (-1 -> 0) = - 

T 

Wi(l -> -1) = ^(0 -)> -1) = - 

T 

Wi (0 -> 1) = -H) = - 



exp(-/3Aj) 


2cosh((35) + exp(- 


-PAi) 


exp(—(35) 




2cosh((35) + exp(- 




exp((35) 




2cosh(f35) + exp(- 





where 5 = JJ2{j) Sj + H. At this point one can notice that Wi(Si id Si new ) 
does not depend on the value Sioid, we can write Wi(Si id — > Si new ) = 
Wi(Sinew), then the master equation becomes: 

—P(S 1: S N ; t) = —J2i J2si i d ^s inew w%(Si new) P (Si, Si old, ■■■Sn', t) 
at ^ 

+ Hi Wi(Si old) J^Si oM ^Si new P(Si, ■■■, Si ne w, ...Sn', t) . 

On the other hand, the sum of probabilities is normalized to one so that 
by multiplying both sides of Eq.Q by S p and taking the average, one obtains, 

rf t (s p ) = - (s p) + ( [ P (Ai) - ; g ^ i 11 u _ J Ai ) . ( i 0) 



2 sinh (B 


J Sj + H 


) 


2cosh(fi 




S, + H 


) + exp(-l3A t ) 
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Finally, after integration over the distribution of -P(Aj) and making use of 
mean field approximation, the kinetic equation of the magnetization becomes, 

„d 2sinh( m+h ^ ) ) , . 2sinh( m+h ™ s ® ) nv 

\l—m = —m + V -r 77^ -, h (1 — P) rr — ttx A ' 

di 2cosh( m+hc T os{i) ) + exp(-|) V ' 2cosh( m+ % osi0 ) + 1 

where £ = ut, m = (S), T = (/3zJ) _1 , d = A/zJ, h = H /zJ . In these 
equations the variable Q was defined as the ratio between the external field 
frequency uj and the frequency of spin flipping (/ = 1/r), i.e., fi = wr = 
u/f. Here we consider a cooperatively interacting many-body system, driven 
by an oscillating external perturbation, an oscillating magnetic field so that 
the thermodynamic response of the system, the magnetization, will then also 
oscillate with necessary modifications in its form [31]. Moreover, the time 
dependence of magnetization can be one of two types according to whether 
they have or do not have the property: 

m{i + tt) = -m(f) . (12) 

A solution that satisfies Eq. ffl2l is called symmetric solution; it corresponds 
to a paramagnetic (P) phase. In this solution, the magnetization m(£) oscil- 
lates around the zero value and is delayed with respect to the external field. 
Solutions of the second type, which do not satisfy Eq. (fl2l . are called non- 
symmetric solutions; they correspond to a ferromagnetic (F) phase. In this 
case, the magnetization does not follow the external magnetic field any more, 
but, instead, oscillates around a nonzero value. Eq. ffTTj) is solved numerically 
by using fourth order Runge-Kutta method for fixed values of T, d,p, and Q. 
Throughout this study we have fixed Q = 2tt, J = 1 and z = 4 for a given set of 
parameters and initial values. The results are presented in Figs.l(a)-(c). Here, 
we can see three different solutions: F,P and coexistence of F and P (F+P). 
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In Fig. 1(a), only the symmetric solution is always obtained, and, hence, we 
have a paramagnetic (P) solution; but, in Fig.l(b), only the nonsymmetric 
solutions are found, and we, therefore, have a ferromagnetic (F) solution. One 
can observe from these figures that these solutions do not depend on the initial 
values. On the other hand, in Fig. 1(c), both the F and P phases exist in the 
system this case corresponds to the coexistence solution (F + P). As can be 
seen in Fig. 1(c) explicitly, the solutions depend on the initial values. 

3 DYNAMIC PHASE TRANSITION POINTS AND PHASE DI- 
AGRAMS 

In order to obtain the dynamic phase boundaries between three phases or 
regions in that are given Figs.l(a)-(c), one should calculate the DPT points. 
The DPT points are obtained by investigating the behavior of the average 
magnetization in a period as a function of the reduced temperature. 



Here m(£) is a stable and periodical function. In general our solution sta- 
bilizes after 6000 periods. In this manner, £ can take any value after this 
transient. In the high field and high temperature region time dependent stag- 
gered magnetization follows the reduced external magnetic field within a single 
period which corresponds to vanishing time average of the dynamical order 
parameter (paramagnetic phase). Whereas, at low field values the magneti- 
zation can not fully switch sign in a single period and the time average of 
the magnetization in a period is non zero and consequently ordered or ferro- 




$0+2tt 



/ 



(13) 
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magnetic phase arises. Fig. 2 represents the reduced temperature dependence 
of the average magnetization (M) for various values of magnetic field ampli- 
tude h and crystal field concentration (p) while d = —0.25. In these figures 
arrows denote the transition temperatures. In Fig.2(a), we give the case for 
h = 0.70, p = 0.50. In this case, the system represents re-entrance with two 
sequential first order phase transitions which take place at T tl and T t2 . While 
Fig. 2(b) exhibits the reduced temperature dependence of the dynamical or- 
der parameter for h = 0.70, p = 0.75. For these values of the parameters, 
BC model with random single ion anisotropy undergoes a first and a second 
phase transition sequentially. In Fig. 2(c) we give an example of second or- 
der phase transition from ordered to disordered phase for h = 0.4, p = 0.75. 
Eventually, Fig. 2(c) illustrates an first order phase transition which occurs for 
/i = 0.75,p = 0.75. 

Fig.3(a) illustrates the thermal variations of M for various values of crystal 
field concentration (p) for vanishing external field. The number accompany- 
ing each curve illustrates the value of p. Here the crystal field has negative 
sign (d = —0.3) since then the critical temperature increases with increasing 
diluteness for fixed h. 

On the other hand, it is well known fact that in the static limit (u = 0.0) 
the dynamic transition disappears and the phase boundary in the h — T plane 
collapses to a line with h = and ending at T = T c , the static transition 
temperature of the unperturbed system [31]. Fig.3(b) shows the thermal vari- 
ations of M and for several values of static h while d = —0.5. In addition, 
Fig. 3(c) gives the behaviors of M and as a function of static h for d = —0.5 
and several values of T. One can see from these figures that the system does 
not undergo any phase transitions for static h. Consequently, we can conclude 
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that the oscillating external magnetic field induces the phase transitions. 

We can now focus on the phase diagrams of the system. In Figs. 3-5 we rep- 
resent the calculated mean field dynamical phase diagrams in the (T,h) plane 
which exhibit the effect of the randomness in the crystal field on the five dif- 
ferent phase diagram topologies of the pure kinetic BC model [30]. First in 
Fig.3(a), we give the dynamical phase diagram of the two-dimensional kinetic 
BC model with bimodal crystal field distribution for d = 0.25. The number 
accompanying each curve denotes the value of the crystal field concentration 
(p). The outermost curve corresponds to the pure BC model with no quenched 
randomness (p = 0).As crystal field quenched randomness is introduced with 
decreasing values of p, ordered phases and first-order phase transitions recede. 
This result is consistent with the RG theory predictions given in Ref . [T7| . Fi- 
nally, we should stress that similar phase diagrams were also obtained in the 
kinetic of the mixed spin-| and spin-1 Ising ferromagnetic system [29] as well 
as the kinetic spin-| Ising model [32]. The reason that the phase diagram is 
similar to the one obtained for the kinetic spin-^ Ising model is due to the 
competition between J, d and h. For positive crystal field values, the Hamilto- 
nian of the spin-1 model gives similar results to the Hamiltonian of the spin-| 
Ising model. 

It has been given in Ref. [30] that pure kinetic BC model has four different 
phase diagram topologies for negative d, which depend on d values. Now let 
us discuss the effect of randomness in the single ion isotropy on these phase 
diagrams: 

(1) For —0.0104 > d > —0.4654, the dynamical phase diagram topology 
of the pure kinetic BC model is similar to positive d case but only differs in 
that for very low T and h values, one more P+F coexistence region also exists. 
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The boundary between this F+P region and the F phase is the first-order line 
(see Fig. 7(b) in Ref. [30J. In BC model, described by the Hamiltonian given in 
Eq. ([T]), negative crystal field interaction (d = A/zJ < 0) favors the annealed 
vacancies, namely the nonmagnetic states Si — 0. Fig. 3(b) exhibits this fact: 
with increasing concentration of negative single ion isotropy (d = —0.25) the 
ordered phase recedes and the tricritical temperature moves to lower temper- 
atures. Whereas, the coexistence region (F+P) in the low temperature and 
field region disappears with increasing vacancy in the single ion anisotropy 
(for p < 0.95). 

(2) For —0.4654 < d < —0.5543, the system exhibits two dynamic tricriti- 
cal points (DTCP). One of them occurs in similar places in the phase diagrams 
for d = 0.25 and d = —0.25 , whereas, the other DTCP occurs in the low h 
region. In addition, the first-order phase transition lines exist at the low re- 
duced temperatures, and h values separate not only the P+F region from the 
F phase, but also from the P phase. When we introduce quench disorder in 
the crystal field we found that this topology changes drastically with varying 
crystal field concentration (p). Our calculations has revealed that there are 
four different phase diagram topologies which depend on p values: 

(2. a) Type 1 (p > 0.97): the system has the same phase diagram with the 
pure case (see Fig. 3(c)). 

(2.b) Type 2 (0.97 > p > 0.85): The DTCP in the low temperature and 
field region disappears. Whereas, the second order line intersects the h = 
axis. Moreover, the P+F phase recedes and an ordered phase appears in the 
neighborhood of h = 0, T = (see Fig.3(d)). 

(2.c) Type 3 (0.85 > p > 0.31): with increasing single ion isotropy va- 
cancy, the coexistence region in the low H and low T region moves to higher 
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field values (see Fig. 3(e)). The diagram contains one DTCP first-order phase 
transition lines merge and signals the change from a first- to a second-order 
phase transitions. 

(2.d) Type 4 (0.31 > p > 0.0): Now the crystal field interaction is rather 
diluted. As one can see from Fig. 3(f) that the F+P phase coexistence totally 
disappears and the phase diagram is similar to the —0.0104 > d > —0.4654 
case. 

(3) For —0.5543 < d < —0.9891, pure system exhibits an exotic phase diagram 
which contains, three different the F+P regions at low reduced temperatures 
in addition to the P and F phases and two dynamical tricritial points. On 
the other hand, if one introduces disorder in the crystal field there are four 
different phase diagram topologies which depend on the concentration of the 
crystal field (see Figs. 5 (a)- (d). ): 

(3. a) Type 1 (1.0 > p > 0.8): Although the dynamical phase diagram of 
the random single ion-anisotropy BC model has similar topology with the pure 
kinetic spin-1 BC model, one can observe that (see Fig. 7(d) of Ref.[30] ) the 
ordered phase moves to lower temperatures and the coexistence region in the 
low temperature and field region shrinks with raising amount of the crystal 
field randomness and finally disappears at p=0.79. 

(3.b) Type 2 (0.80 > p > 0.32): For this interval of the crystal field con- 
centration, the ferromagnetic phase expands in the expense of the F+P and 
P phases and due to this effect the first order transition line which takes place 
in the high field low temperature regime disappears and the system has only 
one DTCP. 

(3.c) Type 3 (0.32 > p > 0.1): The phase boundary that separates F+P 
coexistence phases and F phase turns out to be second order as consequence 
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the dynamical TCP moves to zero temperature . Whereas, the boundary be- 
tween F+P and P phases remains as first order. In addition to the DTCP, 
system exhibits a dynamical critical end point (DCEP) which is shown Fig. 
(5.c) 

(3.d) Type 4 (0.1 > p > 0.0): The behavior differs from Type 3 in the 
sense that for very low T and h values, P+F coexistence region disappears 
and system represents no DCEP. For this interval of the concentration value, 
the system has only one DTCP which exists at low temperature and high 
magnetic field. 

(4) For —0.9891 > d, the topology of the phase diagram is dramatically dif- 
ferent from the other intervals of the single ion anisotrpy amplitude : it does 
not include a P+F phase coexistence region at low temperature and low mag- 
netic field for pure spin-1 BC model. In Figs. 6(a) to (d) we illustrate the four 
different types of behavior depending on the the p: 

4(a) Type 1(1 > p > 0.994): For this interval of the crystal field ampli- 
tude, the system exhibits two DTCP's and the topology of the phase diagram 
is very similar to the pure case. 

4(b) Type 2 (0.993 >p> 0.14): The phase boundary that separates F+P 
coexistence phases and F phase turns to be second order line with decreasing 
concentration value of the randomness single-ion anisotropy. As a result of 
this, the undermost DTCP turns out to be a DCEP. 

4(c) Type 3 (0.14 > p > 0): Finally, for d=-l and p=0, the system exists 
one DTCP at low temperature and high magnetic field. The point where the 
two boundary lines merge. Also this behavior is similar to kinetic mixed Ising 
[29J and kinetic spin- 1/2 Ising model [32J. 
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4 SUMMARY AND CONCLUSIONS 



Within the mean field approach, we have analyzed stationary states of 
the spin-1 Blume-Capel model with a random crystal field Aj under a time- 
dependent oscillating external magnetic field. The time evolution of the system 
is described by a stochastic dynamics of the Glauber type. We have studied the 
time dependence of the magnetization and the behavior of the dynamical or- 
der parameter as a function of reduced temperature for reduced magnetic field 
and different possibility (p) of the crystal field. We have also analyzed thermal 
variations and temperature dependence of M for various values of crystal field 
concentration (p) and for different static reduced magnetic field, respectively. 
Moreover, the behavior of M as function of static reduced magnetic field (h) 
for various values of the reduced temperature have been examined. 

The dynamic phase transition (DPT) points are found and the phase di- 
agrams are constructed in the reduced magnetic field and temperature plane. 
We have found that the behavior of the system strongly depends on the val- 
ues of random crystal field or random single-ion anisotropy. For all (p) and 
positive values of reduced crystal field (d) the system behaves as the standard 
kinetic Ising model |32j, and also kinetic mixed Ising spin- (1/2,1) model [29]. 
We have observed that there exist F+P coexistence and first order dynamical 
phase transitions in the low temperature and high field regime. Whereas, the 
dynamical phase transitions turns out to be second order with increasing tem- 
perature and decreasing magnetic field. Consequently, it shows that the system 
exhibits dynamical tricritical point (DTCP). As we have mentioned in detail 
in the previous section, the introduction of random ion-anisotropy in kinetic 
spin-1 BC model produces an effect which suppress the F+P coexistence region 
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in the low temperature low field and high field low temperature regimes. We 
should stress that this result in according with the previous results obtained 
by Renormalization-Group theory [T7j. Finally we should point out that some 
of the first-order phase lines and also the dynamic multicritical points (DTCP 
and DCEP) are very likely artifacts of the mean-field approach. The reason 
of this artifact can be stated as follows: for field amplitudes less than the 
coercive field and temperatures lower than the static ferromagnetic - param- 
agnetic transition temperature, the time dependent magnetization represents 
a nonsymmetric stationary solution even zero frequency limit. Meanwhile in 
the absence of the fluctuations, the system is trapped in one well of the free 
energy and cannot go to other one [33]. On the other hand, this mean- field 
dynamic study reveals that the random single ion anisotropy spin 1 Blume- 
Capel Model represents interesting dynamic phase diagram topologies. Since 
then, we hope that this work can stimulate further studies on kinetic features 
of kinetic random single ion anisotropy spin-1 Model systems theoretically and 
experimentally. 
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Fig. 1. Time variance of the magnetization (m(£)) while p=0.9 : (a) Corresponding 
to a paramagnetic phase (P) for d=-0.25, h=0.5, and T=0.7; (b) Exhibiting a ferro- 
magnetic phase (F) for d=-0.25, h=0.25, and T=0.5; (c) Representing a coexistence 
region (F+P) d=-0.25, h=0.75, and T =0.1. 
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(c) (d) 

Fig. 2. Dynamical order parameter as a function of reduced temperature. T c and 
Tt indicate second and first order phase transition temperatures respectively, (a) 
The system under goes two successive first order phase transitions, there exists 
re-entrance, (b) Two successive phase transitions: the first one is a first-order and 
the second one a continuous phase transition and there is re-entrance, (c) The 
system under goes a second order phase transition, (d) The system shows a first 
order phase transition. 
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Fig. 3. (a) Thermal variations of M for various values of crystal field concentration 
(p) for vanishing external field. The number accompanying each curve illustrates 
the value of p. (b) Temperature dependence of M for several values of static external 
field amplitudes (h) while p=0.5. The number accompanying each curve denotes the 
value of h. (c) The behavior of M as function of static h for d=-0.5. The number 
accompanying each curve denotes the value of the reduced temperature (T). 
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Fig. 4. Dynamic phase diagrams of the Blume-Capel model with crystal field ran- 
domness in the (T,h) plane for various values of the single ion anisotropy concen- 
tration (p). Dotted and solid lines represent the first-order and second-order phase 
transitions, respectively, (a) d=0.25, the number accompanying each curve denotes 
the value of p. (b) d=-0.25 the number accompanying each curve illustrates the 
value of p. (c) d=-0.525 and p=0.94. (d) d=-0.525 and p=0.85. (e) d=-0.525 and 
p=0.75. (f) d=-0.525 and p=0.25. 
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Fig. 5. Dynamic phase diagrams of the Blumc-Capcl model with crystal field ran- 
domness in the (T,h) plane for various values of the single ion anisotropy con- 
centration (p) while d=-0.625. Dotted and solid lines represent the first-order and 
second-order phase transitions, respectively, (a) p=0.95, (b) p=0.50, (c) p=0.25, (d) 
p=0.0. 
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Fig. 6. Dynamic phase diagrams of the Blume-Capel model with crystal field ran- 
domness in the (T,h) plane for different values of the single ion anisotropy con- 
centration (p) while d=-1.0. Dotted and solid lines represent the first-order and 
second-order phase transitions, respectively, (a) p=1.0, (b) p=0.75, (c) p=0.5, (d) 
p=0.0. 
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